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Abstract 

Within the first 20 minutes of the evolution of the hot, dense, early Universe, astrophysically 
interesting abundances of deuterium, helium-3, helium-4, and lithium-7 were synthesized 
by the cosmic nuclear reactor The primordial abundances of these light nuclides produced 
during Big Bang Nucleosynthesis (BBN) are sensitive to the universal density of baryons and 
to the early-Universe expansion rate which at early epochs is governed by the energy density 
in relativistic particles ("radiation") such as photons and neutrinos. Some 380 kyr later, 
when the cosmic background radiation (CBR) radiation was freed from the embrace of the 
ionized plasma of protons and electrons, the spectrum of temperature fluctuations imprinted 
on the CBR also depended on the baryon and radiation densities. The comparison between 
the constraints imposed by BBN and those from the CBR reveals a remarkably consistent 
picture of the Universe at two widely separated epochs in its evolution. Combining these 
two probes leads to new and tighter constraints on the baryon density at present, on possible 
new physics beyond the standard model of particle physics, as well as identifying some 
challenges to astronomy and astrophysics. In this review the current status of BBN will be 
presented along with the associated estimates of the baryon density and of the energy density 
in radiation. 



1.1 Introduction 

The present Universe is observed to be expanding and filled with radiation (the 
2.7 K cosmic background radiation; CBR) as well as with "ordinary matter" (baryons), 
"dark matter," and "dark energy." As a consequence, the early Universe must have been 
hot and dense. Sufficiently early in its evolution, the universal energy density would have 
been dominated by relativistic particles ("radiation dominated"). During its early evolution 
the Universe passed through a brief epoch when it functioned as a cosmic nuclear reactor, 
synthesizing the lightest nuclides: D, ^He, "^He, and ^Li. These relics from the distant past 
provide a unique window on the early evolution of the Universe, as well as being valuable 
probes of the standard models of cosmology and particle physics. Comparing the predicted 
primordial abundances with those inferred from observational data tests the standard models 
and may uncover clues to their modifications and/or to extensions beyond them. It is clear 
that Big Bang Nucleosynthesis (BBN), one of the pillars of modern cosmology, has a crucial 
role to play as the study of the evolution of the Universe enters a new, data-rich era. 



1 



G. Steigman 



As with all science, cosmology depends on the interplay between theoretical ideas and 
observational data. As new and better data become available, models may need to be refined, 
revised, or even replaced. A consequence of this is that any review such as this one is merely 
a signpost along the road to a better understanding of our Universe. While details of the 
current "standard" model, along with some of its more popular variants to be discussed 
here, may need to be revised or rejected in the future, the underlying physics to be described 
here can provide a useful framework and context for understanding those changes. Any 
quantitative conclusions to be reached today will surely need to be modified in the light 
of new data. This review is, then, a status report on the standard model, highlighting its 
successes as well as exposing the current challenges it faces. While we may rejoice in the 
consistency of the standard model, there is still much work, theoretical and observational, to 
be done. 

1.2 An Overview of BBN 

To set a context for the confrontation of theoretical predictions with observational 
data it is useful to review the physics and cosmology of the early evolution of the Uni- 
verse, touching on the specifics relevant for the synthesis of the fight nucUdes during the 
first ~ 20 minutes. In this section is presented an overview of this evolution along with 
the predicted primordial abundances, first in the standard model and then for two exam- 
ples of nonstandard models which involve variations on the early-Universe expansion rate 
(Steigman, Schramm, & Guim 1977) or asymmetries between the number of neutrinos and 
antineutrinos (e.g., Kang & Steigman 1992, and references therein). 

1.2.1 Early Evolution 

Discussion of BBN can begin when the Universe is a few tenths of a second old and 

the temperature is a few MeV. At such an early epoch the energy density is dominated by the 
relativistic (R) particles present, and the Universe is said to be "radiation-dominated." For 
sufficiently early times, when the temperature is a few times higher than the electron rest- 
mass energy, these are photons, pairs, and, for the standard model of particle physics, 
three flavors of left-handed (i.e., one helicity state) neutrinos (and their right-handed an- 
tineutrinos). 
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PR = P7 + Pe + 3p,.= -^p^, (1.1) 

where is the energy density in CBR photons (which, today, have redshifted to become 
the CBR photons at a temperature of 2.7 K). 

In standard BBN (SBBN) it is assumed that the neutrinos are fully decoupled prior to 
annihilation and do not share in the energy transferred from the annihilating pairs to 
the CBR photons. In this approximation, in the post-e^ annihilation Universe, the photons 
are hotter than the neutrinos by a factor T-y/Tv = {\\ /4)'/^, and the relativistic energy density 
is 

PR = P7 + 3p^ = 1.68p^. (1.2) 

During these radiation-dominated epochs the age (?) and the energy density are related by 

^^^PrI^ = 1, so that once the particle content (pr) is specified, the age of the Universe is 
known (as a function of the CBR temperature Zy). In the standard model. 
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Pm-e^ annihilation : t = 0.738 MeV^ s, (1.3) 

Post-e^ annihilation : t = 1.32 MeV^ s. (1.4) 

Also present at these early times are neutrons and protons, albeit in trace amounts com- 
pared to the relativistic particles. The relative abundance of neutrons and protons is deter- 
mined by the charged-current weak interactions. 

p + e~ < > n + ff, n + e'*' < > P + i^e, n < > p + e~+i>i,. (1-5) 

As time goes by and the Universe expands and cools, the lighter protons are favored over 
the heavier neutrons and the neutron-to-proton ratio decreases, initially following the equi- 
librium form (n/p)eq oc exp(-Am/r), where Am = 1.29 MeV is the neutron-proton mass 
difference. As the temperature drops the two-body collisions in Equation 1 1 . 51 become too 
slow to maintain equilibrium and the neutron-to-proton ratio, while continuing to decrease, 
begins to deviate from (exceeds) this equilibrium value. For later reference, we note that if 
there is an asymmetry between the numbers of Ve and ("neutrino degeneracy"), described 
by a chemical potential /i^ (such that for > there are more ly^ than D^), then the equi- 
librium neutron-to-proton ratio is modified to (n/p) cx exp(-Am/T - jj^e/T). In place of 
the neutrino chemical potential, it is convenient to introduce the dimensionless degeneracy 
parameter = He/T, which is invariant as the Universe expands. 

Prior to annihilation, at T « 0.8 MeV when the Universe is ~ 1 second old, the two- 
body reactions regulating n/p become slow compared to the universal expansion rate and 
this ratio "freezes in," although, in reality, it continues to decrease, albeit more slowly than 
would be the case for equilibrium. Later, when the Universe is several hundred seconds 
old, a time comparable to the neutron lifetime (t„ = 885.7 ± 0.8 s), the n/p ratio resumes 
falling exponentially: n/p oc exp(-f/r„). Since there are several billion CBR photons for 
every nucleon (baryon), the abundances of any complex nuclei are entirely negligible at 
these early times. 

Notice that since the n/p ratio depends on the competition between the weak interaction 
rates and the early-Universe expansion rate (as well as on a possible neutrino asymmetry), 
any deviations from the standard model (e.g., pR pR + px or ^ 0) will change the relative 
numbers of neutrons and protons available for building more complex nuclides. 

1.2.2 Building the Elements 

At the same time that neutrons and protons are interconverting, they are also col- 
liding among themselves to create deuterons: n + p < — > D + 7. However, at early times, 
when the density and average energy of the CBR photons are very high, the newly formed 
deuterons find themselves bathed in a background of high-energy gamma rays capable of 
photodissociating them. Since there are more than a billion photons for every nucleon in 
the Universe, before the deuteron can capture a neutron or a proton to begin building the 
heavier nuclides, the deuteron is photodissociated. This bottleneck to BBN persists until 
the temperature drops sufficiently so that there are too few photons energetic enough to pho- 
todissociate the deuterons before they can capture nucleons to launch BBN. This occurs after 
annihilation, when the Universe is a few minutes old and the temperature has dropped 
below 80 keV (0.08 MeV). 

Once BBN begins in earnest, neutrons and protons quickly combine to form D, ^H, ''He, 
and "^He. Here, at "^He, there is a different kind of bottleneck. There are no stable mass-5 
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nuclides. To jump this gap requires "*He reactions with D or -^H or ^He, all of which are 
positively charged. The Coulomb repulsion among these colliding nuclei suppresses the 
reaction rates, ensuring that virtually all of the neutrons available for BBN are incorporated 
in ^He (the most tightly bound of the light nuclides), and also that the abundances of the 
heavier nuclides are severely depressed below that of "^He (and even of D and ^He). Recall 
that is unstable, decaying to ^He. The few reactions that manage to bridge the mass-5 gap 
lead mainly to mass-7 (^Li or ^Be, which, later, when the Universe has cooled further, will 
capture an electron and decay to ^Li); the abundance of ^Li is 1 to 2 orders of magnitude 
below that of the more tightly bound ^Li. Finally, there is another gap at mass-8. This 
absence of any stable mass-8 nuclei ensures there will be no astrophysically interesting 
production of any heavier nuclides. 

The primordial nuclear reactor is short-Uved, quickly encountering an energy crisis. Be- 
cause of the falling temperature and the Coulomb barriers, nuclear reactions cease rather 
abruptly as the temperature drops below ~ 30 keV, when the Universe is ~ 20 minutes old. 
This results in "nuclear freeze-out," since no already existing nuclides are destroyed (except 
for those that are unstable and decay) and no new nuclides are created. In ^ 1000 seconds 
BBN has run its course. 

1.2.3 The SBBN -predicted Abundances 

The primordial abundances of D, ■'He, and ^Li(^Be) are rate limited, depending 
sensitively on the competition between the nuclear reaction rates (proportional to the nu- 
cleon density) and the universal expansion rate. As a result, these nuclides are all potential 
baryometers. As the Universe expands, the nucleon density decreases so it is useful to com- 
pare it to that of the CBR photons: 77 = WN/n^. Since this ratio turns out to be very small, it 
is convenient to introduce 

r?io = 10'"(nN/«7) = 274r!b/z2 , (1 .6) 

where fib is the ratio of the present values of the baryon and critical densities and h is the 
present value of the Hubble parameter in units of 100 km s~' Mpc"' . As the Universe 
evolves (post-e^ annihilation) this ratio is accurately preserved so that 77 at the time of BBN 
should be equal to its value today. Testing this relation over ten orders of magnitude in 
redshift, over a timespan of some 10 billion years, can provide a confirmation of, or pose a 
challenge to the standard model. 

In contrast to the other light nuclides, the primordial abundance of "^He (mass fraction 
Y) is relatively insensitive to the baryon density, but since virtually all neutrons available 
at BBN are incorporated in '*He, Y does depend on the competition between the weak in- 
teraction rates (largely fixed by the accurately measured neutron lifetime) and the universal 
expansion rate. The higher the nucleon density, the earlier can the D bottleneck be breached. 
Since at early times there are more neutrons (as a fraction of the nucleons), more '*He will be 
synthesized. This latter effect is responsible for a very slow (logarithmic) increase in Y with 
77. Given the standard model relation between time and temperature and the measured nu- 
clear and weak cross sections and decay rates, the evolution of the light-nuclide abundances 
may be calculated and the relic, primordial abundances predicted as a function of the one 
free parameter, the nucleon density or 77. These predictions for SBBN are shown in Figure 

o 
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Fig. 1.1. The SBBN-predicted primordial abundances of D, ^He, and ^Li (by number with 
respect to hydrogen), and the *He mass fraction Y as a function of the nucleon abundance 
7710. The widths of the bands reflect the theoretical uncertainties. 

Not shown on Figure FTTTI are the relic abundances of ^Li, ^Be, "^B, and "B; for the same 
range in ?/, all of them lie offscale, in the range 10"^"- 10"'^. The results shown here are 
from the BBN code developed and refined over the years by my colleagues at The Ohio State 
University (OSU). They are in excellent agreement with the published results of the Chicago 
group (Buries, Nollett, & Turner 2001). Notice that the abundances appear in Figure fLTI as 
bands. These reflect the theoretical uncertainties in the predicted abundances. For the OSU 
code the errors in D/H and ^He/H are at the ^ 8% level, while they are much larger, ~ 12%, 
for ^Li. Buries et al. (2001), in a reanalysis of the relevant published cross sections, have 
reduced the theoretical errors by roughly a factor of 3 for D and ^He and a factor of 2 for ^Li. 
The reader may not notice the band shown for "^He, since the uncertainty in Y, dominated by 
the very small uncertainty in the neutron lifetime, is at only the ^ 0.2% level (ay ~ 0.0005). 
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Based on the discussion above it is easy to understand the trends shown in Figure [TTTI D 
and ^He are burned to "^He. The higher the nucleon density, the faster this occurs, leaving 
behind fewer nuclei of D or ■'He. The very slight increase of Y with 77 is largely due to 
BBN starting earlier at higher nucleon density (more complete burning of D, ^H, and ^He 
to ^He) and higher neutron-to-proton ratio (more neutrons, more "^He). The behavior of ^Li 
is more interesting. At relatively low values of ryio 1S 3, mass-7 is largely synthesized as 
^Li [by ^H(Q;,7)^Li reactions], which is easily destroyed in collisions with protons. So, as 77 
increases at low values, destruction is faster and ^Li/H decreases. In contrast, at relatively 
high values of 7710 ^ 3, mass-7 is largely synthesized as ^Be [via ''He(a,7)^Be reactions], 
which is more tightly bound than ^Li and, therefore, harder to destroy. As rj increases at 
high values, the abundance of ^Be increases. Later in the evolution of the Universe, when it 
is cooler and neutral atoms begin to form, ^Be will capture an electron and /3-decay to ^Li. 

1.2.4 Nonstandard BBN 

The predictions of the primordial abundance of ^He depend sensitively on the early 
expansion rate (the Hubble parameter H) and on the amount — if any — of a z/^ - P,, asym- 
metry (the I'e chemical potential fi^ or the neutrino degeneracy parameter ^g). In contrast 
to ^He, the BBN-predicted abundances of D, ^He and ^Li are determined by the competi- 
tion between the various two-body production/destruction rates and the universal expansion 
rate. As a result, the D, ^He, and ^Li abundances are sensitive to the post-e='= annihilation 
expansion rate, while that of "'He depends on both the pre- and post-e^ annihilation expan- 
sion rates; the former determines the "freeze-in" and the latter modulates the importance 
of /3-decay (see, e.g., Kneller & Steigman 2003). Also, the primordial abundances of D, 
^He, and ^Li, while not entirely insensitive to neutrino degeneracy, are much less affected 
by a nonzero (e.g., Kang & Steigman 1992). Each of these nonstandard cases will be 
considered below. Note that the abundances of at least two different relic nuclei are needed 
to break the degeneracy between the baryon density and a possible nonstandard expansion 
rate resulting from new physics or cosmology, and/or a neutrino asymmetry. 

1.2.4.1 Additional Relativistic Energy Density 

The most straightforward variation of SBBN is to consider the effect of a non- 
standard expansion rate H' ^ H. To quantify the deviation from the standard model it is 
convenient to introduce the "expansion rate factor" (or speedup/slowdown factor) S, where 

S = H'/H = t/t'. (1.7) 

Such a nonstandard expansion rate might result from the presence of "extra" energy con- 
tributed by new, light (relativistic at BBN) particles "X". These might, but need not, be 
additional flavors of active or sterile neutrinos. For X particles that are decoupled, in the 
sense that they do not share in the energy released by annihilation, it is convenient to 
account for the extra contribution to the standard-model energy density by normalizing it to 
that of an "equivalent" neutrino flavor (Steigman et al. 1977), 

px = /:^upu='^-^^N^p^. (1.8) 

For SBBN, AA^i. = (A^^ = 3 + /S.N,j) and for each such additional "neutrino-like" particle 
(i.e., any two-component fermion), if Tx = T^,, then AA^i, = l;ifX should be a scalar, ANi, = 
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All . However, it may well be that the X have decoupled even earlier in the evolution of the 
Universe and have failed to profit from the heating when various other particle-antiparticle 
pairs annihilated (or unstable particles decayed). In this case, the contribution to /S.Ny from 
each such particle will be < 1 (< 4/7). Henceforth we drop the X subscript. Note that, 
in principle, we are considering any term in the energy density that scales like "radiation" 
(i.e., decreases with the expansion of the Universe as the fourth power of the scale factor). In 
this sense, the modification to the usual Friedman equation due to higher dimensional effects, 
as in the Randall-Sundrum model (Randall & Sundrum 1999a,b; see also Cline, Grojean, & 
Servant 1999; Binetruy et al. 2000; Bratt et al. 2002), may be included as well. The interest 
in this latter case is that it permits the possibility of an apparent negative contribution to 
the radiation density {iSNi, < 0; 5 < 1). For such a modification to the energy density, the 
pre-e^ annihilation energy density in Eauation ll.ll is changed to 

(Pr),..= + (1-9) 

Since any extra energy density (JS.Nu > 0) speeds up the expansion of the Universe (S > 
1), the right-hand side of the time-temperature relation in Equation 11.31 is smaller by the 
square root of the factor in parentheses in Equation ll.9l 

43 

In the post-e^ annihilation Universe the extra energy density is diluted by the heating of 
the photons, so that 



Spre = (t/t')pre = (1 + -7^^)'^' = ( 1 + 0. 163 AA?,)'/^ (I.IQ) 



and 



iPR)pos, = 1.68(1+0.1 35AN,)p^ (1.11) 



Spos, = {t/t')posr = (l+0.135AA^,)'/2. (1.12) 



While the abundances of D, ''He, and ^Li are most sensitive to the baryon density (77), the 
^He mass fraction (Y) provides the best probe of the expansion rate. This is illustrated in 



Figure fL2l where. in the AA^^ - 7710 plane, are shown isoabundance contours for D/H and Yp 
(the isoabundance curves for ^He/H and for ^Li/H, omitted for clarity, are similar in behavior 
to that of D/H). The trends illustrated in Figure [l~2l are easy to understand in the context 
of the discussion above. The higher the baryon density (ryio), the faster primordial D is 
destroyed, so the relic abundance of D is anticorrelated with 7710. But, the faster the Universe 
expands (AA^^ > 0), the less time is available for D destruction, so D/H is positively, albeit 
weakly, correlated with AA^i,. In contrast to D (and to ^He and ^Li), since the incorporation 
of all available neutrons into "^He is not limited by the nuclear reaction rates, the '^He mass 
fraction is relatively insensitive to the baryon density, but it is very sensitive to both the 
pre- and post-e^ annihilation expansion rates (which control the neutron-to-proton ratio). 
The faster the Universe expands, the more neutrons are available for '^He. The very slow 
increase of Yp with 7710 is a reflection of the fact that for a higher baryon density, BBN begins 
earUer, when there are more neutrons. As a result of these complementary correlations, the 
pair of primordial abundances = 10^(D/H)p and Yp, the ^He mass fraction, provide 
observational constraints on both the baryon density (77) and on the universal expansion rate 
factor S (or on AA^^,) when the Universe was some 20 minutes old. Comparing these to 
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Fig. 1.2. Isoabundance curves for D and "^He in the AA^^ - 7710 plane. The solid curves are 
for '*He (from top to bottom: Y = 0.25, 0.24, 0.23). The dotted curves are for D (from left 
to right; = 10^(D/H) = 3.0, 2.5, 2.0). The data point with error bars corresponds to = 
2.6±0.4 and Yp = 0.238±0.005; see the text for discussion of these abundances. 



similar constraints from when the Universe was some 380 Kyr old, provided by the WMAP 
observations of the CBR polarization and the spectrum of temperature fluctuations, provides 
a test of the consistency of the standard models of cosmology and of particle physics and 
further constrains the allowed range of the present-Universe baryon density (e.g., Barger et 
al. 2003a,b; Crotty, Lesgourgues, & Pastor 2003; Hannestad 2003; Pierpaoli 2003). 

1.2.4.2 Neutrino Degeneracy 

The baryon-to-photon ratio provides a dimensionless measure of the universal baryon 
asymmetry, which is very small (77 < 10""^). By charge neutrality the asymmetry in the 
charged leptons must also be of this order. However, there are no observational constraints, 
save those to be discussed here (see Kang & Steigman 1992; Kneller et al. 2002, and further 
references therein), on the magnitude of any asymmetry among the neutral leptons (neu- 
trinos). A relatively small asymmetry between electron type neutrinos and antineutrinos 
(Cc ^ 10"^) can have a significant impact on the early-Universe ratio of neutrons to protons, 
thereby affecting the yields of the light nuclides formed during BBN. The strongest effect 
is on the BBN "^He abundance, which is neutron limited. For S,e > 0, there is an excess of 
neutrinos (vg) over antineutrinos (Dg), and the two-body reactions regulating the neutron-to- 
proton ratio (Ea. ll.5> drive down the neutron abundance; the reverse is true for < 0. The 
effect of a nonzero asymmetry on the relic abundances of the other light nuclides is much 
weaker. This is illustrated in Figure [OJ which shows the D and '^He isoabundance curves 
in the - 7710 plane. The nearly horizontal "^He curves reflect the weak dependence of Yp on 



8 



G. Steigman 



0.1 



0) 
JJLT 




-0.05 ' — ' — ' — ' — ' — ^ — ' — ' — — ' — ' — — ' — ' — ' — ^ — ' — ' — ' — ' 
4 5 6 7 8 

Tlio 



Fig. 1.3. Isoabundance curves for D and ^He in the -?7io plane. The solid curves are for 
"^He (from top to bottom: Yp = 0.23, 0.24, 0.25). The dotted curves are for D (from left to 
right: y-Q = 10^(D/H) = 3.0, 2.5, 2.0.) The data point with error bars corresponds to = 
2.6±0.4 and Yp = 0.238±0.005; see the text for discussion of these abundances. 

the baryon density, along with its significant dependence on the neutrino asymmetry. In con- 
trast, the nearly vertical D curves reveal the strong dependence of on the baryon density 
and its weak dependence on any neutrino asymmetry (^He/H and ^Li/H behave similarly: 
strongly dependent on ry, weakly dependent on ^p). This complementarity between y^ and 
Yp permits the pair [rj^^e] to be determined once the primordial abundances of D and ^He 
are inferred from the appropriate observational data. 

1.3 Primordial Abundances 

It is clear from Figures [l~n - ll.3l that tests of the consistency of SBBN, along with 
constraints on any new physics, will be data-driven. While D (and/or ^He and/or ^Li) largely 
constrain the baryon density and "^He plays a similar role for AA^^ and/or for ^p, there is an 
interplay among 7710, /SN^, and which is quite sensitive to the adopted abundances. For 
example, a lower primordial D/H increases the BBN-inferred value of 7710, leading to a 
higher predicted primordial ^He mass fraction. If the primordial ^He mass fraction derived 
from the data is "low," then a low upper bound on AA^^ (or a nonzero lower bound on ^g) 
will be inferred. It is therefore crucial to avoid biasing any conclusions by underestimating 
the present uncertainties in the primordial abundances derived from the observational data. 

The four light nuclides of interest, D, ^He, '*He, and ^Li follow very different evolutionary 
paths in the post-BBN Universe. In addition, the observations leading to their abundance 
determinations are also very different. Neutral D is observed in absorption in the UV; singly 
ionized ^He is observed in emission in Galactic H 11 regions; both singly and doubly ionized 
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"^He are observed in emission via recombinations in extragalactic H ii regions; ^Li is observed 
in absorption in the atmospheres of very metal-poor halo stars. The different histories and 
observational strategies provide some insurance that systematic errors affecting the inferred 
primordial abundances of any one of the light nuclides are unlikely to distort the inferred 
abundances of the others. 

1.3.1 Deuterium 

The post-BBN evolution of D is straightforward. As gas is incorporated into stars 
the very loosely bound deuteron is burned to ^He (and beyond). Any D that passes through a 
star is destroyed. Furthermore, there are no astrophysical sites where D can be produced in 
an abundance anywhere near that observed (Epstein, Lattimer, & Schramm 1976). As a re- 
sult, as the Universe evolves and gas is cycled through generations of stars, deuterium is only 
destroyed. Therefore, observations of the deuterium abundance anywhere, anytime, provide 
lower bounds on its primordial abundance. Furthermore, if D can be observed in "young" 
systems, in the sense of very little stellar processing, the observed abundance should be 
very close to the primordial value. Thus, while there are extensive data on deuterium in the 
solar system and the local interstellar medium of the Galaxy, it is the handful of observa- 
tions of deuterium absorption in high-redshift, low-metallicity QSO absorption-line systems 
(QSOALS), which are potentially the most valuable. At sufficiently high redshifts and low 
metallicities, the primordial abundance of deuterium should reveal itself as a "deuterium 
plateau." 

Inferring the primordial D abundance from the QSOALS has not been without its dif- 
ficulties, with some abundance claims having been withdrawn or revised. Presently there 
are ~ half a dozen QSOALS with reasonably firm deuterium detections (Buries & Tytler 
1998a,b; D'Odorico, Dessauges-Zavadsky, & Molaro 2001; O'Meara et al. 2001; Pettini 
& Bowen 2002; Kirkman et al. 2003). However, there is significant dispersion among the 
derived abundances, and the data fail to reveal the anticipated deuterium plateau (Fig. 11.41 
- 11.61 see also Steigman 2003). Furthermore, subsequent observations of the D'Odorico 
et al. (2001) QSOALS by Levshakov et al. (2002) revealed a more complex velocity struc- 
ture and led to a revised — and uncertain — deuterium abundance. This sensitivity to poorly 
constrained velocity structure in the absorbers is also exposed in the analyses of published 
QSOALS data by Levshakov and collaborators (Levshakov, Kegel, & Takahara 1998a,b, 
1999), which lead to consistent, but somewhat higher, deuterium abundances than those 
inferred from "standard" data reduction analyses. 

Indeed, the absorption spectra of D i and H i are identical, except for a wavelength/velocity 
offset resulting from the heavier reduced mass of the deuterium atom. An H i "interloper," 
a low-column density cloud shifted by ^ 81 km s"' with respect to the main absorbing 
cloud, would masquerade as Di . If this is not accounted for, a D/H ratio which is too 
high would be inferred. Since there are more low-column density absorbers than those 
with high H i column densities, absorption-line systems with somewhat lower H i column 
density (e.g., Lyman-lknit systems) are more susceptible to this contamination than are the 
higher H i column density absorbers (e.g., damped Lya absorbers). However, for the damped 
Lya absorbers, an accurate determination of the H i column density requires an accurate 
placement of the continuum, which could be compromised by interlopers. This might lead 
to an overestimate of the Hi column density and a concomitant underestimate of D/H (J. 
Linsky, private communication). As will be seen, there is the possibility that each of these 
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Fig. 1.4. The deuterium abundance, D/H, versus metallicity, "X"(usually, X = Si), from 
observations (as of early 2003) of QSOALS (filled circles). Also shown for comparison 
are the D abundances for the local ISM (filled square) and the solar system ("Sun"; filled 
triangle). 



effects may have contaminated the current data. Indeed, complex velocity structure in the 
D'Odorico et al. (2001) absorber (see Levshakov et al. 2002) renders it of less value in 
constraining primordial deuterium, and it will not be included in the estimates presented 
here. 

In Figure [l~4l are shown the extant data (circa June 2003) for D/H as a function of metal- 
Hcity from the work of Buries & Tytler (1998a,b), O'Meara et al. (2001), Pettini & Bowen 
(2002), and Kirkman et al. (2003). Also shown for comparison are the local interstellar 
medium (ISM) D/H (Linsky & Wood 2000) and that for the presolar nebula as inferred from 
solar system data (Geiss & Gloeckler 1998). 

On the basis of our discussion of the post-BBN evolution of D/H, a "deuterium plateau" at 
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Fig. 1.5. The deuterium abundance, D/H, vers us the Hi column density in the absorbers, 
N(H I ), for the same QSOALS as in FigurelTH 



low metallicity was expected. If, indeed, one is present, it is hidden by the dispersion in the 
current data. Given the possibility that interlopers may affect both the D i and the H i column 
density determinations, it is interesting to plot D/H as a function of N(Hi ). This is shown 
in Figure [O] where there is some (limited) evidence that D/H is higher in the Lyman-limit 
systems than in the damped Lya absorbers. 

To decide how to utilize this confusing data it may be of value to consider the observa- 
tions chronologically. Of the set chosen here. Buries & Tytler (1998a,b) studied the first 
two lines of sight. For PKS 1937-1009 they derived yo = 105(D/H) = 3.25 ±0.3 (Buries 
& Tytler 1998a), while for Q1009+299 they found = 3-9S^di (Buries & Tytler 1998b). 
These two determinations are in excellent agreement with each other (x^ = 1.0), leading 
to a mean abundance (yu) = 3.37 ±0.27. Next, O'Meara et al. (2001) added the line of 
sight to HS 0105+1619, finding a considerably lower abundance yu = 2.54 ±0.23. In- 
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Fig. 1.6. As in Figure \L4\ The dashed Unes represent the ±la band calculated from the 
mean and its dispersion ((D/H)p = 2.6 ±0.4 x 10"^; see the text). 



deed, while the weighted mean for these three hnes of sight is (yo) = 2.88, the has 
ballooned to 6.4 (for two degrees of freedom). Absent any evidence that one or more of 
these abundances is in error, O'Meara et al. adopt the mean, and, for the error in the mean, 
they take the dispersion about the mean (0.72) divided by the square root of the number 
of data points: (yo) = 2.88 ±0.42. One year later Pettini & Bowen (2002) published their 
//ST data on the line of sight toward Q2206-199, finding a surprisingly low value of yr, = 
1 .65 ± 0.35. Including this determination reduces the mean to (yo) = 2.63, but the disper- 
sion in grows to 1.00 and = 16.3 for three degrees of freedom. Clearly, either one 
or more of these determinations is in error, or the variation among the high-redshift, low- 
metallicity deuterium abundances is larger than anticipated from our understanding of its 
evolution (Jedamzik & Fuller 1997). Using the mean and its dispersion (to fix the error), as 
of the time of the Carnegie Symposium, the best estimate for the primordial D abundance 
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was (yu) = 2.63 ±0.50. Shortly thereafter, in early 2003, the data of Kirkman et al. (2003) 
appeared for the line of sight toward Q 1243+3047. For this hne of sight they find yu = 
2.42;!;q 25. This abundance lies between the lowest and the higher previous values, reducing 
the overall dispersion to 0.88, while hardly changing the mean fromyo = 2.63 to 2.60. While 
the total is still enormous, increasing slightly to 16.6, the reduced decreases from 5.4 
to 4.2. This is still far too large, suggesting that one or more of these determinations may be 
contaminated, or that there may actually be real variations in D/H at high redshifts and low 
metallicities. Notice (see Fig. ll.5> that the largest D/H estimates are from the two absorbers 
with the lowest H i column densities (Lyman-limit systems), where interlopers might con- 
tribute to the inferred D i column densities, while the lowest abundances are from the higher 
H I column density (damped Lya ) absorbers, where interlopers might affect the wings of the 
H I lines used to fix the H i column densities. Absent any further data supporting, or refuting, 
these possibilities, there is no a priori reason to reject any of these determinations. 

To utilize the current data, the weighted mean D abundances for these five lines of sight 
and the dispersion are used to infer the abundance of primordial deuterium (and its uncer- 
tainty) adopted in this review: yu = 2.6 ± 0.4. Note that, given the large dispersion, two- 
decimal place accuracy seems to be wishful thinking at present. For this reason, in quoting 
the primordial D abundance inferred from the observational data 1 have purposely chosen to 
quote values to only one decimal place. This choice is consistent too with the ^ 3%-8% 
theoretical uncertainty (at fixed i]) in the BBN-predicted abundance. In Figure [l~6l are shown 
the data, along with the corresponding la band. It is worth remarking that using the same 
data Kirkman et al. (2003) derive a slightly higher mean D abundance: yo = 2.1 A. The rea- 
son for the difference is that they first find the mean of log(yD) and then use it to compute 
the mean D abundance (jd = lO^'"^*-''^))). 



1.3.2 Helium-3 

The post-BBN evolution of ^He is considerably more complex and model depen- 
dent than that of D. Interstellar ^He incorporated into stars is burned to ^He (and beyond) in 
the hotter interiors, but preserved in the cooler, outer layers. Furthermore, while hydrogen 
burning in cooler, low-mass stars is a net producer of ^He (Iben 1967; Rood 1972; Dear- 
born, Schramm, & Steigman 1986; Vassiliadis & Wood 1993; Dearborn, Steigman, & Tosi 
1996) it is unclear how much of this newly synthesized -^He is returned to the interstellar 
medium and how much of it is consumed in post-main sequence evolution (e.g., Sackmann 
& Boothroyd 1999a,b). Indeed, it is clear that when the data (Geiss & Gloeckler 1998; Rood 
et al. 1998; Bania, Rood, & Balser 2002) are compared to a large variety of chemical evolu- 
tion models (Rood, Steigman, & Tinsley 1976; Dearborn et al. 1996; GalH et al. 1997; Palla 
et al. 2000; Chiappini, Renda, & Matteucci 2002), agreement is only possible for a very 
delicate balance between net production and net destruction of -^He. For a recent review of 
the current status of ^He evolution, see Romano et al. (2003). Given this state of affairs it is 
not possible to utilize ■'He as a baryometer, but it may perhaps be used to provide a consis- 
tency check. To this end, the abundance inferred by Bania et al. (2002) from an H ii region 
in the outer Galaxy, where post-BBN evolution might have been minimal, is adopted here: 
y3 =105(^He/H) =1.1 ±0.2. 
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1.3.3 Helium-4 

Helium-4 is the second most abundant nuclide in the Universe after hydrogen. In 
post-BBN evolution gas cychng though stars has its hydrogen burned to helium, increasing 
the ^He abundance above its primordial value. As with deuterium, a ^He "plateau" is ex- 
pected at sufficiently low metallicity. Although "^He is observed in the Sun and in Galactic 
H II regions, the crucial data for inferring its primordial abundance is from observations of 
the helium and hydrogen emission (recombination) lines from low-metallicity, extragalac- 
tic Hii regions. The present inventory of such regions studied for their helium content is 
approaching of order 100. Thus, it is not surprising that even with modest observational 
errors for any individual Hii region, the statistical uncertainty in the inferred primordial 
abundance may be quite small. In this situation, care must be taken with hitherto ignored or 
unaccounted for corrections and systematic errors or biases. 

In Figure [l~7l is shown a compilation of the data used by Olive & Steigman (1995) and 
Olive, Skillman, & Steigman ( 1 997), along with the independent data set obtained by Izotov, 
Thuan, & Lipovetsky (1997) and Izotov & Thuan (1998). To track the evolution of the "^He 
mass fraction, Y is plotted versus the Hii region oxygen abundance. These Hii regions are 
all metal poor, ranging from ^ \/2 down to '--^ 1 /40 of solar (for a solar oxygen abundance 
of O/H « 5 X 10"^; Allende-Prieto, Lambert, & Asplund 2001). A key feature of FigureO 
is that for sufficiently low metallicity the Y versus O/H relation approaches a "^He plateau! 
Since Y increases with metallicity, the relic abundance can either be bounded from above 
by the lowest metallicity regions, or the Y versus O/H relation may be extrapolated to zero 
metallicity. The extrapolation is quite small, so that whether the former or the latter approach 
is adopted the difference in the inferred primordial abundance is small: | AY| ^ 0.001. 

While the data shown in Figure [LTl reveal a well-defined primordial abundance for '^He, 
the scale hides the very small statistical errors as well as the tension between the two 
groups' helium abundances. Olive & Steigman (1995) and Olive et al. (1997) find Yp = 
0.234±0.003,but Izotov et al. (1997) and Izotov & Thuan (1998) derive Yp = 0.244 ±0.002. 
Although it is difficult to account for all of the difference, much of it is traceable to the dif- 
ferent ways the two groups correct for the contribution to the emission lines from coUisional 
excitation of neutral helium and also to Izotov and collaborators rejecting some hehum emis- 
sion lines a posteriori when they yield "too low" an abundance. Furthermore, for either data 
set, there are additional corrections for temperature, for temperature and density fluctuations, 
and for ionization, which when applied can change the inferred primordial ^He abundance 
by more than the quoted statistical errors (see, e.g., Steigman, Viegas, & Gruenwald 1997; 
Viegas, Gruenwald, & Steigman 2000; Gruenwald, Steigman, & Viegas 2002; Peimbert, 
Peimbert & Luridiana 2002; Sauer & Jedamzik 2002). 

For example, Peimbert et al. (2002) recently reanalyzed the data from four of the Izotov 
& Thuan (1998) Hii regions, employing their own Hii region temperatures and accounting 
for temperature fluctuations. Peimbert et al. (2002) derive systematically lower helium 
abundances, as shown in Figure [l~8l From this very limited sample Peimbert et al. suggest 
that the Izotov & Thuan (1998) estimate for the primordial "^He mass fraction might have to 
be reduced by as much as ^ 0.007. Peimbert et al. go further, combining their redetermined 
helium abundances for these four H ii regions with an accurate determination of Y in a more 
metal-rich Hii region (Peimbert, Peimbert, & Ruiz 2000). Although these five data points 
are consistent with zero slope in the Y - O/H relation, leading to a primordial abundance 
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Fig. 1.7. The "^He mass fraction, Y, inferred from observations of low-metallicity, extra- 
galactic H II regions versus the oxygen abundance derived from the same data. (Figure cour- 
tesy of K. A. Olive.) 

Yp = 0.240±0.001, this extremely small data set is also consistent with AY« 40(O/H), 
leading to a smaller primordial estimate of Yp w 0.237. 

It seems clear that until new data address the unresolved systematic errors afflicting the 
derivation of the primordial helium abundance, the true errors must be much larger than the 
statistical uncertainties. In an attempt to account for this, here I follow Olive, Steigman, & 
Walker (2000) and adopt a compromise mean value along with a larger uncertainty: Yp = 
0.238 ±0.005. 

1.3.4 Lithium-7 

Lithium-7 is fragile, burning in stars at a relatively low temperature. As a result, 
the majority of any interstellar ^Li cycled through stars is destroyed. For the same reason, it 
is difficult for stars to create new ^Li and/or to return any newly synthesized 'Li to the ISM 
before it is destroyed by nuclear burning. In addition to synthesis in stars, the intermediate- 
mass nuclides ^Li, 'Li, ^Be, '°B, and "B can be synthesized via cosmic ray nucleosynthesis, 
either by alpha-alpha fusion reactions, or by spallation reactions (nuclear breakup) in colli- 
sions between protons and alpha particles and CNO nuclei. In the early Galaxy, when the 
metallicity is low, the post-BBN production of lithium is expected to be subdominant to that 
from BBN abundance. As the data in Figure [L9l i'eveal. only relatively late in the evolution 
of the Galaxy does the lithium abundance increase. The data also confirm the anticipated 
"Spite plateau" (Spite & Spite 1982), the absence of a significant slope in the Li/H versus 
[Fe/H] relation at low metallicity due to the dominance of BBN-produced 'Li. The plateau 
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Fig. 1.8. The Peimbert et al. (2002) reanalysis of the "^He abundance data for four of the 
Izotov & Thuan (1998) Hii regions. The open circles are the Izotov & Thuan (1998) abun- 
dances, while the filled circles are from Peimbert et al. (2002). 



is a clear signal of the primordial lithium abundance. Notice, also, the enormous spread 
among the lithium abundances at higher metallicity. This range in Li/H likely results from 
the destruction/dilution of lithium on the surfaces of the observed stars while they are on the 
main sequence and/or lithium destruction during their pre-main sequence evolution, imply- 
ing that it is the upper envelope of the Li/H versus [Fe/H] relation that preserves the history 
of Galactic lithium evolution. Note, also, that at low metallicity the dispersion is much nar- 
rower, suggesting that corrections for depletion/dilution are (may be) much smaller for the 
Population 11 stars. 

As with the other relic nuclides, the dominant uncertainties in estimating the primordial 
abundance of ^Li are not statistical, but systematic. The lithium observed in the atmospheres 
of cool, metal-poor. Population 11 halo stars is most relevant for determining the BBN ^Li 
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Fig. 1.9. A compilation of the lithium abundance data as a function of metallicity from 
stellar observations (courtesy of V. V. Smith). e(Li) = 10'^(Li/H), and [Fe/H] is the usual 
logarithmic metallicity relative to solar Note the "Spite plateau" in Li/H for [Fe/H] ^ - 2. 



abundance. Uncertainties in the lithium equivalent width measurements, in the temperature 
scales for the cool Population 11 stars, and in their model atmospheres dominate the overall 
error budget. For example, Ryan et al. (2000), using the Ryan, Norris, & Beers (1999) 
data, infer [Li]p = 12-l-log(Li/H) = 2.1, while Bonifacio & Molaro (1997) and Bonifacio, 
Molaro, & Pasquini (1997) derive [Li]p = 2.2, and Thorburn (1994) finds [Li]p = 2.3. From 
recent observations of stars in a metal-poor globular cluster, Bonifacio et al. (2002) derive 
[Li]p = 2.34 ± 0.056. As may be seen from Figure [L9l the indication from the preliminary 
data assembled by V. V. Smith (private communication) favors a Spite plateau at [Li]p « 2.2. 

In addition to these intrinsic uncertainties, there are others associated with stellar structure 
and evolution. The metal-poor halo stars that define the primordial lithium plateau are very 
old. As a result, they have had time to disturb the prestellar lithium that could survive in their 
cooler, outer layers. Mixing of these outer layers with the hotter interior where lithium has 
been (can be) destroyed will dilute or deplete the surface lithium abundance. Pinsonneault 
et al. (1999, 2002) have shown that rotational mixing may decrease the surface abundance 
of lithium in these Population 11 stars by 0.1 - 0.3 dex while still maintaining the rather 
narrow dispersion among the plateau abundances (see also Chaboyer et al. 1992; Theado 
& Vauclair 2001; Salaris & Weiss 2002). Pinsonneault et al. (2002) adopted for a baseline 
(Spite plateau) estimate [Li] = 2.2 ±0.1, while for an overall depletion factor 0.2 ±0.1 
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dex was chosen. Adding these contributions to the log of the primordial lithium abundance 
linearly, an estimate [Li]p = 2.4 ±0.2 was derived. In the comparison between theory and 
observation below, I will adopt the Ryan et al. (2000) estimate [Li]p = 2.1 ±0.1, but 1 will 
also consider the implications of the Pinsonneault et al. (2002) value. 



1.4 Confrontation of Theory with Data 

Having reviewed the basic physics and cosmological evolution underlying BBN 
and summarized the observational data leading to a set of adopted primordial abundances, 
the predictions may now be confronted with the data. There are several possible approaches 
that might be adopted. The following option is chosen here. First, concentrating on the 
predictions of SBBN, deuterium will be used as the baryometer of choice to fix the baryon- 
to-photon ratio ij. This value and its uncertainty are then used to "predict" the ^He, "^He, and 
^Li abundances, which are compared to those adopted above. This comparison can provide a 
test of the consistency of SBBN as well as identify those points of "tension" between theory 
and observation. This confrontation is carried further to consider the two extensions beyond 
the standard model [S i 1 (AA?^ t^O); ^ 0]. 



1.4.1 Testing the Standard Model 

For SBBN, the baryon density corresponding to the D abundance adopted here 
(yo = 2.6 ±0.4) is 7/10 = 6.1^0 5, corresponding to fib = 0.022^Q Qg2- This is in outstanding 
agreement with the estimate of Spergel et al. (2003), based largely on the new CBR (WMAP) 
data (Bennett et al. 2003): fib = 0.0224 ± 0.0009. For the baryon density determined by D, 
the SBBN-predicted abundance of ^He is 3^3 = 1.0 ±0.1, which is to be compared to the 
outer-Galaxy abundance of 3/3 = 1.1 ±0.1, which is suggested by Bania et al. (2002) to be 
nearly primordial. Again, the agreement is excellent. 

The tension between the data and SBBN arises with "^He. Given the very slow variation 
of Yp with J], along with the very high accuracy of the SBBN-predicted abundance, the pri- 
mordial abundance is tightly constrained: Ysbbn = 0.248 ± 0.001 . This should be compared 
with our adopted estimate of Y = 0.238 ± 0.005 (Olive et al. 2000). Agreement is only at 
the ^ 5% level. This tension is shown in Figure [TTTOI This apparent challenge to SBBN is 
also an opportunity. As already noted, while the "*He abundance is insensitive to the baryon 
density, it is very sensitive to new physics (i.e., nonstandard universal expansion rate and/or 
neutrino degeneracy). 

There is tension, too, when comparing the SBBN-predicted abundance of ^Li with the 
(very uncertain) primordial abundance inferred from the data. For SBBN the expected abun- 
dance is [Li]p = 2.65^0 [J. This is to be compared with the various estimates above that 
suggested [Li]p w 2.2 ±0.1. In Figure fLTTl is shown the analog of Figure [TTTOI for lithium 
and deuterium. Depending on the assessment of the uncertainty in the primordial abundance 
inferred from the observational data, the conflict with SBBN may or may not be serious. In 
contrast to "^He, ^Li is more similar to D (and to ^He) in that its BBN-predicted abundance is 
relatively insensitive to new physics. As a result, this tension, if it persists, could be a signal 
of interesting new astrophysics (e.g., have the halo stars depleted or diluted their surface 
lithium?). 
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Fig. 1.10. The SBBN-predicted relation between the primordial abundances of D and ^He 
(mass fraction) is shown by the band, whose thickness represents the uncertainties in the 
predicted abundances. Also shown by the point and error bars are the adopted primordial 
abundances of D and "^He (see the text). 



1.4.2 Nonstandard Expansion Rate: Sj^ I ( ANi, y^Q) 

The excellent agreement between the SBBN-predicted baryon density inferred from 
the primordial-D abundance and that derived from the CBR and large scale structure (Spergel 
et al. 2003), and also the agreement between predicted and observed D and ''He suggest that 
the tension with ^He, if not observational or astrophysical in origin, may be a sign of new 
physics. As noted earlier, Yp is sensitive to the early-Universe expansion rate (while D, ''He, 
and 'Li are less so). A faster expansion (S > 1, ANi, > 0) leads to a higher predicted pri- 
mordial abundance of ^He, and vice versa for S < I (ANi, < 0). In Figure [l~T2l is shown the 
same Yp versus yu band as for SBBN in Figure lLTOI along with the corresponding bands for 
the nonstandard cases of a faster expansion {ANi, = 4) and a slower expansion (AA^i, = 2). 
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Fig. 1.11. The SBBN-predicted relation between the primordial abundances of D and ^Li is 
shown by the band, whose thickness reflects the uncertainties in the predicted abundances. 
The data points are for the primordial abundance of D adopted here along with the Ryan et 
al. (2000) Li abundance (filled circle) and the Pinsonneault et al. (2002) Li abundance (filled 
triangle). 



It can be seen that the data "prefer" a slower than standard early-Universe expansion rate. 
If both rj and ANi, are allowed to be free, it is possible (not surprisingly) to accommodate 
the adopted primordial abundances of D and '^He (see Fig. \1.2i . Given the similar effects 
of AN^ T^'O on the BBN-predicted D, ^He, and ^Li abundances, while it is possible to main- 
tain the good agreement (from SBBN) for ^He, the tension between ^Li and D cannot be 
relieved. In Fi sure [TTSl are shown the 1-, 2-, and 3-a BBN contours in the rj - AN,^ plane 
derived from the adopted values of yu and Yp. Although the best-fit point is at AA^,^ = -0.7 
(and 7710 = 5.7), it is clear that SBBN (Ni, = 3) is acceptable. 

The CBR temperature anisotropy spectrum and polarization are also sensitive to the early- 
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Fig. 1.12. As in Figure lTTOI foi- = 2, 3, 4, which coiTespond to S = 0.915, 1, 1.078. 



Universe expansion rate (see, e.g., Barger et al. 2003a, and references therein). There is 
excellent overlap between the 77 - AA^i, confidence contours from BBN as shown in Figure 
11.131 and from the CBR (Barger et al. 2003a). In Figure 11.141 are shown the confidence 
contours in the rj - ANi, plane for a joint BBN - CBR fit (Barger et al. 2003a). Again, while 
the best fit value for ANi, is negative (driven largely by the adopted value for Yp ), AA^^ = 
is quite acceptable. 

1.4.3 Neutrino Asymmetry (i^^ j^Q) 

The tension between D and '*He can also be relieved by nonstandard neutrino 
physics (see Fig. II. 3> . Although the asymmetry (difference between the numbers of par- 
ticles and antiparticles) in charged leptons, tied to that in the baryons by charge neutrality of 
the Universe, must be very small, the neutrino asymmetry is unconstrained observationally. 
Of relevance to BBN is the asymmetry between the electron neutrinos and the electron an- 
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tineutrinos (^g), which regulates the pre-BBN neutron-to-proton ratio through the reactions 
in Eauation ll.51 In Figure lTTTSl are shown the 1- and 2-a contours in the 77-^^ plane for BBN 
(for A^i, = 3) and the adopted abundances of D and ^He. As seen before for AA^^ j^O, while 
a fit to the data can be achieved for j^O, the data are not inconsistent with = 0. Further- 
more, as is shown in Figure lTTTSl BBN constrains the allowed range for neutrino asymmetry 
to be very small. For further implications for neutrino physics and for a discussion of the 
case where both ANi, and are free to differ from zero, see Barger et al. (2003b). 

1.5 Summary and Conclusions 

Given the standard models of cosmology and particle physics, SBBN predicts the 
primordial abundances of D, ^He, ^He, and ^Li, which may be compared with the observa- 
tional data. Of the light nuclides, deuterium is the baryometer of choice, while '^He is an ex- 
cellent chronometer. The universal density of baryons inferred from SBBN and the adopted 
primordial D abundance is in excellent (exact!) agreement with that derived from non-BBN, 
mainly CBR data (Spergel et al. 2003): ?7io(SBBN) = 6.10!^;^; r7io(CBR) = 6.14 ±0.25. 
For this baryon density, the predicted primordial abundance of ^He is also in excellent 
agreement with the (very uncertain) value inferred from observations of an outer-Galaxy 
Hn region (Bania et al. 2002). In contrast, the SBBN-predicted mass fraction of "^He for 
the concordant baryon density is Yp = 0.248 ±0.001, while that inferred from observations 
of recombination lines in metal-poor, extragalactic H 11 regions is lower (Olive et al. 2000): 
Yp*' = 0.238 ±0.005. Since the uncertainties in the observationally inferred primordial value 
are likely dominated by systematics, this ~ 2a difference may not be cause for (much) con- 
cern. Finally, there appears to be a more serious issue concerning the predicted and observed 
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Fig. 1 . 14. The 1 - and 2-a contours in the 77 - AA^^, plane for the joint BBN - CBR (WMAP) 
fit (Barger et al. 2003a). 



lithium abundances. While the predicted abundance is [Li]p « 2.6 ±0.1, current observa- 
tions of metal-poor halo stars suggest a considerably smaller value « 2.2 ± 0. 1 . 

It has been seen that the tension between D and ^He (or between the baryon density and 
^He) can be relieved by either of two variations of the standard model (slower than standard 
early expansion rate; nonzero chemical potential for the electron neutrino). However, in 
neither of these cases does the BBN-predicted ^Li abundance move any closer to that inferred 
from the observations. 

In the current, data-rich era of cosmological research, BBN continues to play an important 
role. The spectacular agreement in the baryon density inferred from processes occurring at 
widely separated epochs confirms the general features of the standard models of cosmology 
and particle physics. The tensions with "^He and ^Li provide challenges, and opportuni- 
ties, to cosmology, to astrophysics, and to particle physics. To outline these challenges and 
opportunities, let us consider each of the light nuclides in turn. 

For deuterium the agreement between SBBN and non-BBN determinations is perfect. 
This may be surprising given the unexpectedly large dispersion among the handful of extant 
D abundance determinations at high redshifts and low metallicities. Here, the challenge is to 
observers and theorists. Clearly more data are called for Perhaps new data will reduce the 



24 



G. Steigman 
0.2 I- 




.Q -| I I I I I I I I I I I I I I I I I I I I I I I I I 

4 5 6 7 8 9 

TllO 

Fig. 1.15. The 1-, 2-, and 3-cr contours in the plane for BBN (A^^ = 3) and the adopted 
D and "^He abundances (Barger et al. 2003b). 

dispersion. In that case it can be anticipated that the SBBN-predicted baryon density will 
approach the accuracy of that currently available from non-BBN data. On the other hand, 
newer data may support the dispersion, suggesting unexpectedly large variations in the D 
abundance at evolutionary times earlier than expected (Jedamzik & Fuller 1997). Perhaps 
there is more to be learned about early chemical evolution. 

From studies of ^He in Galactic Hii regions (Balser et al. 1997; Bania et al. 2002) it 
appears that in the course of Galactic chemical evolution there has been a very delicate 
balance between post-BBN production and destruction. If either had dominated, a gradient 
of the ^He abundance with galactocentric distance should have been seen in the data (see 
Romano et al. 2003, and references therein). So far, none is. Clearly, more data and a 
better understanding of the lower mass stars, which should dominate the production and 
destruction of ^He, would be of value. 

The very precise value of the baryon density inferred either from D and SBBN or from 
non-BBN data, coupled with the very weak dependence of the SBBN abundance of ^He 
on the baryon density, leads to a very precise prediction of its primordial mass fraction. 
Although there exists a very large data set of '^He abundance determinations, the observa- 
tional situation is confused at present. It seems clear that while new data would be valuable, 
quality is much more important than quantity. Data that can help resolve various correc- 
tions for temperature, for temperature and density fluctuations, for ionization corrections, 
would be of greater value than merely collecting more data that are incapable of addressing 
these issues. Because of the very large data set(s), the statistical uncertainty in the derived 
primordial mass fraction is very small, avp ~ 0.002-0.003, while uncertain systematic cor- 
rections are much larger }^ 0.005. At this point it is systematics, not statistics that dominate 
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the uncertainty in the primordial helium abundance. In this context it is worth considering 
non-emission line observations that might provide an independent abundance determination. 
Just such an alternative, the so-called R-parameter method using globular cluster stars was 
proposed long ago by Iben (1968) and by Iben & Faulkner (1968). It too has many sys- 
tematic uncertainties associated with its application, but they are different from those for 
the emission-line studies. Very recently, Cassisi, Salaris, & Irwin (2003), using new stellar 
models and nuclear reactions rates, along with better data, find Yp = 0.243 ± 0.006. This 
is in much better agreement with the expected value (within ^ la) and should stimulate 
further investigations. 

The apparent conflict between the predicted and observed abundances of ^Li, if not simply 
traceable to the statistical and systematic uncertainties, suggests a gap in our understanding 
of the structure and evolution of the very old, metal-poor, halo stars. It would appear from 
the comparison between the predicted and observed abundances that lithium may have been 
depleted or diluted from the surfaces of these stars by ^ 0.2-0.4 dex. Although a vari- 
ety of mechanisms for depletion/dilution exist, the challenge is to account for such a large 
reduction without at the same time producing a large dispersion around the Spite plateau. 

The wealth of observational data accumulated over the last decade or more have propelled 
the study of cosmology from youth to maturity. BBN has played, and continues to play, a 
central role in this process. There have been many successes, but much remains to be done. 
Whether the resolution of the current challenges are observational or theoretical, the future 
is bright. 
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